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It is shown that small-scale magnetic fields present before recombination induce baryonic density 
inhomogeneities of appreciable magnitude. The presence of such inhomogeneities changes the ion- 
ization history of the Universe, which in turn decreases the angular scale of the Doppler peaks and 
increases Silk damping by photon diffusion. This unique signature could be used to (dis)prove the 



existence of primordial magnetic fields of strength as small as B \ 
microwave background observations. 
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Primordial magnetic fields may well have been gener- 
ated during early cosmic phase transitions, during an in- 
flationary epoch (in case conformal invariance is broken), 
or during an epoch of baryogenesis, among other. In fact, 
it seems unlikely that the early Universe was not magne- 
tized due to the multitude of possibilities. The question 
is rather of which strength such fields would be, and on 
what typical length scales they would reside. Recently, 
it has been claimed that the surprisingly weak flux of 
GeV 7-rays in the direction of three TeV-blazars may be 
best understood by the presence of cosmic magnetic fields 
of relatively weak strength [l[ filling a large fraction of 
space [3] (see Ref. [3J however). Such fields could explain 
why secondary GeV 7-rays, induced by TeV 7-rays pair 
producing on the infrared background, with the produced 
subsequently inverse Compton scattering on the cos- 
mic microwave background radiation (CMBR), would be 
moved out of the light cone due to the curved trajectories 
of the e . 

It would be important to find other observational 
signatures of such putative primordial magnetic fields. 
A prime candidate here are precision observations of 
anisotropies of the CMBR. A larger number of studies 
have been presented, with the majority of studies assum- 
ing substantial fields on 10-100 Mpc scales [J| (an excep- 
tion is Ref. (H). Such fields can, however, realistically 
only be produced during an inflationary scenario, with 
the stringent requirements of breaking conformal invari- 
ance and avoiding backreaction of the created magnetic 
fields on the inflationary process. To be observable, field 
strength of B ~ 10~ 9 Gauss [6| have to be assumed. 
Field of that strength may, however, already potentially 
be ruled out due to likely overproduction of magnetic 
fields in galaxy clusters Q. 

When magnetogenesis happens after inflation, result- 
ing magnetic field spectra are blue, with much more 
power on small scales than on large scales. For dynam- 
ically relaxed magnetic fields a correlation between the 
final present day magnetic field strength B and its cor- 



relation length L may be given Q 
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On larger scales fields are likely falling of with a white 
noise spectrum B ~ (l/L)~ 3 / 2 Q or even steeper 0. 
Magnetic fields on kpc scales are usually not believed to 
change the observable anisotropies in the CMBR since 
that scale would correspond to multipoles of I ~ 10' 
whereas the Planck mission will observe only up to I ~ 
2 — 3 x 10 3 . We will show here that this view is incorrect, 
i.e. magnetic fields on such small scales do change the 
anisotropies in the CMBR on smaller multipoles. 

Shortly before recombination CMBR photons do not 
participate in fluid flows on kpc scales, as the photon 
mean free path is much larger £ 7 ~ Mpc. They do, how- 
ever, strongly affect fluid flows by introducing a high drag 
on moving electrons due to occasional Thomson scatter- 
ings, leaving the plasma on small scales in a highly vis- 
cous state [l(| before recombination. Immediately after 
the decoupling of photons on scale L (i.e. when Z 7 be- 
comes larger than L) the plasma experiences an enor- 
mous ~ 3 x 10~ 5 d ecrease in the speed of sound from 
eg = 1/ \/3(l + R) where R = 3pb / 4p 7 with pb , p 7 
the photon and baryon mass densities, respectively, to 
c s = v /2T/m b [nj]. That is, whereas for all purposes the 
plasma had been incompressible when l 1 < L it becomes 
compressible, at least for sufficiently large magnetic field 
strength, when Z 7 > L. 

Imagine a stochastic magnetic field and negligible ve- 
locities v initially. The evolution of velocities and densi- 
ties are given by the Euler and continuity equations 
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where a ~ l/i 7 (cf. 7]) is the photon drag term. In the 
overdamped, highly viscous state before recombination, 
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only the terms on the RHS of Eq. ([2} are important. Very 
quickly (At ^ 1/a) terminal velocities of v ~ v\ / aL are 
reached. Here va = B / ^/Airp is the Alfven velocity of the 
baryon plasma. For a stochastic field the generated fluid 
flows are necessarily both rotational (i.e. V x v 7^ 0) 
and compressional (i.e. V • v 7^ 0). In fact, one can 
show that the energy dissipation rate of compressional 
modes is a factor (va/c s ) 2 smaller than those for rota- 
tional modes, such that when c s ^> va a larger part of 
magnetic field configuration which induce compressional 
flows could potentially survive. The compressional com- 
ponent leads to the creation of density fluctuations. Us- 
ing Eq. ^ one finds Sp/p(t) ~vt/L~ v 2 A t/aL 2 . These 
density fluctuations become larger with time until cither, 
pressure forces become important in counteracting fur- 
ther compression, or the source magnetic stress term de- 
cays. The former happens when the last term on the 
LHS of Eq. ((2J (cg/L)Sp/p is of the order of the mag- 
netic force term v\jL. That is, density fluctuations may 
not become larger than Sp/p< (va/c s ) 2 . Here Alfven- 
and sound speed shortly before recombination are given 
by 
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It has been shown in Ref. 0] that magnetic fields do de- 
cay even in the viscous photon free-streaming limit appli- 
cable shortly before recombination. Here decay of mag- 
netic energy occurs via the excitation of fluid flows which 
are than converted to heat due to photon drag. Though 
counterintuitive, dissipation is stronger when the drag 
term a becomes weaker. By direct numerical simula- 
tion the linear analysis [lol | and non-linear estimate [l2j ] 
was confirmed that magnetic fields do decay when the 
eddy turnover rate v/L ~ v 2 A /aL 2 equals the Hubble 
rate H ~ 1/t. Entering this into the above expression 
for 5p/p(t) one finds that the average density fluctua- 
tion is not expected to exceed unity by much, even for 
vanishing c s . Putting all this together, we expect 
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for the density fluctuations generated by magnetic fields 
before recombination. 

The cosmological hydrogen recombination process is 
well approximated by the following differential equation 
in time fl3ll 



dn c 
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where n e , n#o, and nn are electron-, neutral hydrogen- 
and total hydrogen- nn = n^a + n e density respectively, 
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FIG. 1: Ionization fraction X e for a Universe with the 
preferred WMAP7-year cosmological parameters citeWMAP 
as a function of redshift for a homogeneous Universe 
(solid) and two Universes with small-scale inhomogeneities 
y/{(5n/n) 2 ) = y/?> (dotted) and 10 (dashed), respec- 
tively. Here, somewhat arbitrarily, the parameters f v = 
(0.567,0.33,0.1) and A* = (0.1,1,6) for the first inho- 
mogeneous model and f v = (0.663,0.33,0.04) and A i = 
(0.1, 1, 168) for the second model have been chosen, where 
fv is the volume fraction filled with density n 1 = (1 + A J )(n) 
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and with a e , f3 e , and A the Case B recombination rate, 
photoionization rate from the 2s level, and the 2s — > Is 
two photon decay rate, respectively. Furthermore hv a 
is the Lyman-a transition energy, T is temperature, and 
K = \^/(8irH) with A Q the Lyman-a wavelength and 
H the Hubble constant. Note that Eq. ([7]) is only for 
illustrative purposes as it neglects the presence of helium. 

A key observation in Eq. ([7]) is that it is non-linear in 
density in the first term on the RHS and in the factor 
C [HI- In an Universe inhomogeneous on scales L < ( 7 
CMBR anisotropies depend on the average electron den- 
sity (n e ). However, due to the non-linearity (n e ) 7^ nj! omo 
where n^ omo is the electron density in a homogeneous 
Universe, irrespective of the fact that the baryon density 
n homo _ e q ua i s th e average baryon density in the in- 
homogeneous Universe. This may be seen in Fig. 1 where 
we computed with help of the public code RECFAST [l6[ 
the ionization fraction X e = (n e ) / (nn) 7^ {X e ) in a inho- 
mogeneous Universe by taking the average of the electron 
densities of three independent regions [171 ] with different 
baryonic densities but with the same average density as 
a homogeneous Universe. It is seen that the drop in X e , 
i.e. recombination, occurs earlier when inhomogeneities 
exist. 

CMBR temperature anisotropies may be calculated at 
linear order by evolving temperature Oq and gravita- 



tional potential "J perturbations of wavevector k across 
the epoch of recombination. When this is done the ob- 
served temperature fluctuations are related to 9o + * = 
(6 + $)V(k) pj] where 

V(k) = /"dij^e-^'et^M 2 (9) 
Jo d? 7 

is due to the imperfect coupling of photons and 
baryons inducing exponential damping of perturbations 
on the photon diffusion scale ^0(77), and g(n)dn — 
dT/dryexp(— r)d?7 is the fraction of photons observed to- 
day (i.e. at present conformal time rjo) which scattered 
last between conformal times rj and r/ + dn, with 17(77) the 
visibility function. Here dr/dr/ = X e nH<JThi with gtk 
the Thomson cross section and a scale factor, such that 
r is the photon optical depth. The damping factor T>(k) 
strongly modifies the undamped 0o+ v I / temperature fluc- 
tuations. The behavior of Oo is given by the solutions of 
the equation of a forced oscillator. It is well known that 
due to well specified initial conditions (i.e. only grow- 
ing modes) Oo(fc) exhibits an oscillatory behavior with 
peaks given by k m r s ~ rmr where m is an integer. These 
peaks are due to perturbations having performed half, 
one, one-and-a-half, ... sonic oscillations where 

rs=J Q c s d V = J z j sm (10) 

is the sound horizon, with ?y rcc (zdrcc) conformal time 
(redshift) at recombination, c s the baryon-photon speed 
of sound, and H the Hubble constant. Correspond- 
ing peaks in the temperature-temperature correlation 
function on angular scale, or equivalently on spherical 
harmonic multipole I, are observed at l m ~ k m (r/ — 
rj rec ). The above gives us most of the ingredients to 
qualitatively understand modifications in the CMBR 
anisotropies from small scale inhomogeneity. 

We modified CAMB pi to compute the CMBR 
anisotropies when magnetic field induced baryon den- 
sity fluctuations are present. The results are shown in 
Fig. [2] showing the anisotropies and their fractional de- 
viations from the best-fit WMAP7 model for the two 
inhomogeneous Universes with X e shown in Fig. [T] It 
is seen that inhomogeneities have two main effects (a) 
they move the Doppler peak locations to higher multi- 
poles and (b) they enhance Silk damping of the high I 
peaks. Both may be understood by inspecting Fig. ([1}. 
In small-scale inhomogeneous Universes high density re- 
gions recombine earlier, making the average ionization 
fraction X e drop significantly earlier, and therefore in- 
creasing the redshift of ionization z rcc . Low density re- 
gions which recombine later are not too important for the 
visibility function since they do not contain too many 
electrons. For example, for y/ ((Sn/n) 2 ) — -s/3 recom- 
bination (the peak of the visibility function) is moved 
from z rec ~ 1078 to z rcc rj 1118 a substantial change 
of 6.6%. In order to induce such a large change of z Tec 
by a change of the baryonic- or matter- densities, rela- 
tive changes of ~ 20%, ~ 10%, respectively are required. 
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FIG. 2: CMBR anisotropies in conventional units (cf. [HI ]) 
(upper panel) as a function of multipole for the best-fit 
WMAP7 year model, and the two inhomogeneous models de- 
scribed in the caption of Fig. [TJ as well as their fractional 
differences (lower panel) to the best-fit WMAP7 year model. 
The same line coding as in Fig. [T] is applied. 

These are far beyond the WMAP7 and baryonic acous- 
tic oscillation (BAO) [2(| error bars on their respective 
values of ~ ±2%, and ~ ±3.3%, respectively. An earlier 
recombination leads to Doppler peaks moving to higher 
I, e.g. for yf ((Sn/n) 2 ) = y3 all lower peaks are moved by 

~ 2.5%. This may be understood since to lowest order 

1/2 

Al ~ Az r ec (see below). 

The second effect, enhanced Silk damping, is some- 
what more surprising since the Silk damping scale is the 
diffusion scale (i.e. c? 7 ~ y/Tyt with Z 7 photon mean free 
path and t time) at recombination. Earlier recombina- 
tion would imply less time for photon diffusion and so 
less Silk damping. However, inspecting again Fig. ([T]) 
one observes that in the inhomogeneous Universes (X e ) 
is smaller by ~ 10% already some time before recombina- 
tion. This is due to earlier helium recombination in the 
high density regions. A smaller electron density implies 
larger Z 7 and therefore larger d 1 . As this latter effect 
dominates, the combined effect is more Silk damping as 
evident from Fig. @. 

Fractional changes in the CMBR anisotropies in in- 
homogeneous Universes compared to homogeneous Uni- 
verses are substantial ~ 10% even for yj ((Sn/n) 2 ) ~ 1, 
particular at high multipoles. On first sight one would 
think that such large changes may be detected or ruled 
out by the Planck mission. To good approximation l m ~ 
k m f]o ~ Vo/rs ~ Vo/Vrec ~ [H ICC (l + z Q )]/[H (l + 2: rec )]. 
Assume for the moment that the Universe only contains 
matter Q,Mh 2 = fldh 2 + flbh 2 and radiation, where Qd, &b 
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are dark matter- and baryon- contributions to the criti- 
cal density today, and h is the Hubble constant in units 
of 100km s _1 Mpc _1 . If one then assumes that the Uni- 
verse is critically closed l m becomes a function of only z rec 

1/2 

(since [11^(1 + z )]/[H (l + z rcc )] ~ 2tec up to calculable 
radiation contributions) independent of VLm or the Hub- 
ble constant. Since z rcc is only logarithmically depen- 
dent on well-know atomic physics, magnetic field induced 
density fluctuations could be detected/ruled out to very 
high precision via the shift of the Doppler peaks. Un- 
fortunately the situation is somewhat more complicated 
as the present Universe is dominated by a cosmological 
constant. In that case 1/Hq ^770 = fi(flMh 2 , ft Ah 2 ) and 
H r0 c = f2{Qhih 2 , Zrec) (using that the radiation density 
is well known), where fi are functions. To break possible 
degeneracies between density inhomogeneities and other 
cosmological parameters fti^h 2 and f^A^ 2 must be known 
accurately. Here SJ^ft 2 can be inferred from the CMBR 
anisotropics and fl\h 2 ) from either supernovae surveys 
or BAO observations of the angular diameter distance. 
Nevertheless, it is not clear if one can achieve the desired 
accuracy. Alternatively, assuming a closed Universe and 
using a precise measurement of the Hubble constant h 
could also lead to a fairly precise prediction of z rcc and 
the Z m 's. A more detailed analysis is beyond the scope of 
this letter and is deferred to a future publication. 

Which field strengths are detectable in case the Planck 
mission combined with other surveys will be able to es- 
tablish (or refute) the existence of small-scale inhomo- 
geneity before recombination [2lJ ? Following Ref. Q pri- 
mordial magnetic fields do decay on the scale implicitly 



given by v/L ~ H , with v ~ v 2 A / La before recombination 
in the viscous regime and v ~ va after in the turbulent 
regime. Assuming that the initial spectrum is given by 
B ~ L~ n l 2 one may deduce that magnetic field strength 
shortly before, B < and after B > recombination (which 
is also the final present day field strength Q) are related 
by B< ~ {a/H ICC ) n /<- 2n+ ^B>. That is, due to the rapid 
disappearance of photon drag during recombination, sub- 
stantial amounts of magnetic field energy density dissi- 
pates right at recombination since a/H rcc « 170. It is 
not clear, but subject to further investigation, if shocks 
resulting from the magnetic stress acceleration and the 
accompanying shock ionization during recombination are 
also of importance. In any case, for a final magnetic field 
strength of B > ~ 10~ n Gauss and a white noise spectral 
index n = 3 one finds _B < ~ 4.7 x 10~ n Gauss which 
comes close to fulfilling va ~ c s such that it should be 
potentially detectable. 

In summary, we have argued that primordial magnetic 
fields induce small-scale baryon inhomogeneity of sub- 
stantial amplitude. Due to non-linearities in the recom- 
bination equations such inhomogeneity induces changes 
in the ionization fraction before recombination. This, 
in turn, influences the anisotropies in the CMBR, such 
that present day primordial magnetic fields of strength 
B ~ 10~ n Gauss could potentially be detected by a com- 
bination of future and present CMBR, BAO, and super- 
novae observations. 
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